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Abstract For the first time a quasi-simultaneous multi- 
wavelength campaign has been performed on an Anoma- 
lous X-ray Pulsar from the radio to the hard X-ray band. 
4U 0142+61 was an INTEGRAL target for 1 Ms in July 
2005. During these observations it was also observed in 
the X-ray band with Swift and RXTE, in the optical 
and NIR with Gemini North and in the radio with the 
WSRT. In this paper we present the source-energy dis- 
tribution. The spectral results obtained in the individual 
wave bands do not connect smoothly; apparently com- 
ponents of different origin contribute to the total spec- 
trum. Remarkable is that the INTEGRAL hard X-ray 
spectrum (power-law index 0.79 ±0.10) is now measured 
up to an energy of ^230 keV with no indication of a spec- 
tral break. Extrapolation of the INTEGRAL power-law 
spectrum to lower energies passes orders of magnitude 
underneath the NIR and optical fluxes, as well as the 
low ~30 fj,3y (2a) upper limit in the radio band. 
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1 Introduction 

Anomalous X-ray Pulsars (AXPs) are young rotati ng iso- 
lated neutron stars (for a review in this volume, see lKaspj 
120061 ). Currently ther e are 8 AXPs known and the re are a 
few more candidates (|Woods fc Thompsonl[2006] ). These 
objects are called anomalous, because their X-ray lumi- 
nosities exceed by far the available total energy released 
by rotational energy loss. The energy output is believed 
to originate from an immense energy reservoir stored in a 
toroidal magnetic field within the Neutron Star. The sur- 
face magnetic fields, inferred from their periods and pe- 
riod derivatives, are of the order oflO 14 — 10 15 G. There- 
fore, AXPs are believed to be magnetars, as originally 
proposed for the Soft Gamma-ray Repeaters (SGRs, see 
Duncan fc Thompsol]|l992trTn~ompson fc Duncanlll995l 
Il996t ). Both AXPs and SGRs are well studied objects in 
the X-ray band for energies below 10 keV. However, little 
was known about their persistent emission in the hard 
X-ray band (>10 keV). In 2004 INTEGRAL discovered 
hard X-ray emissio n from the position of IE 1841-045 
(|Molkov et alJ2004h . lKuiper. Hermsen fc Mende3 (|2004 ) 
showed unambiguously that the hard X-rays originated 
from the AXP by extracting a pulsed hard X-ray signal 
from the source using archival RXTE data. After INTE- 
GRAL discovered hard X-rays from tw o other AXP loca- 
tions , namely from 1R XS J1708-4009 (iRevniytsev et al.l 
I2004D and 4U 0142+61 (den Hartog et al.ll2004f i. lKuiper et all 
(2006}) also showed for these and for a fourth AXP (IE 2259+586) 



pulsed hard X-ray emission using archival RXTE data. 
That means that presently already for 4 of the 8 estab- 
lished AXPs hard X-ray emission has been detected and 
this can now be considered to be a common characteris- 
tic, which is not yet understood. 
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In this paper we focus on the AXP 4U 0142+61. This 
AXP was discovered by the Uhuru X-ray observatory in 
the e arly seventies ( Giacconi et al.l Il972t iForman et al.l 
Il978l) . The spin period of 8.7 s was found by Israel, 
Mereghetti & Stella (1994). They realised that the X- 
ray luminosity is too high to be explained by rotational 
energy loss. Like for the other AXPs, there is no proof for 
a companion, nor for an active (i.e. accreting) disk that 
could explain the high X-ray luminosity of 4U 0142+61. 
T he passive (i.e. non accreting) debris disk discovered 
bv lWang. Chakrabartv k Kapl an ( 2006) does not power 
the X-ray emission. The X-ray luminosities recently mea- 
sured with e.g. XMM-Newton and Chandr a are of the 
order of 10 35 erg cm" 2 s" 1 (2-10 j keV, see iPatel et all 
120031: iGohler. Wilms k Staubertl 2005 )), assumi ng a dis- 
tance of 3.6 kpc (|Durant k van Kerkwiildl2006al) . The X- 
ray spectra (0.5-10 keV) of AXPs are soft and are com- 
monly fitted with a black-body and a power-law model. 
The inclusion of the power-law component is required 
to fit excess photons with energies above ~3 keV. For 
4U 0142+61, the best fit parameters are a black-body 
temperature of kT ~0.4 keV and a power-law photon 
index of r ~3.4. 



5<t upper limit of 0.27 mJy could be extracted from the 
data. 

4U 0142+61 was the first AXP for which an optical 
counterpart was discovered (Hul l eman, van Kerkwijk & 
Kulkarni 2000). iKern k Martml (|2002f l found that the 
optical emission is pulsed for a considera ble fraction. 
iHulleman. van Kerkwiik k Kulkarnil ( 2000T ) showed for 
the first time that it is not possible to understand the 
optical and NIR measured fluxes with respect to the X- 
ray fluxes. There is an optical and NIR excess that can 
not be explained by the Rayleigh- Jeans tail from the X- 
ray black body. Moreover the optical and NIR emissions 
seem to be non ther mal and exhibit more variability than 
seen in the X-rays dHulle man, van Kerkwiik k Kulkarnil 



2004llsrael et al l2004lMorii et al.l2005HDurant k van Kerkwiik! 



4U 0142+61 was detected bv lden Hartog et al.l (|2006ah 
in hard X-rays up to 1 50 keV in 1.6 Ms o f INTEGRAL 
observations (see also iKuiper et al. 2006). The 20 150 
keV flux was measured to be (9.7+0.9) x 10 -11 erg cm -2 s~ 
The total spectrum could be fitted with a power-law 
model with photon index r — 0.73 ± 0.17. They also re- 
visited the Compton Gamma-Ra y Observatory (CQMP- 
TEL) archives (0.75-30 MeV, see lSchonfelder et al.ll2002D 
and determined flux upper limits at the location of the 
AXP. These limits put constraints on the extrapolation 
of the hard X-ray power law. Assuming that the hard 
X-ray flux is stable, a spectral break has to occur in 
the hard X-ray regime below ^750 keV in order for the 
spectrum n ot to be in conflic t with t he C OMPTEL up- 
per li mits. Iden Hartog et ail (2006a) and IKuiper et al.l 
(2006) were not able to measure such a break, but a hint 
was found with a 2.3cr fit improvement when a high- 
energy cutoff was added to the power-law model by den 
Hartog et al. (2006a). The indication for the cutoff was 
found at an energy of 73 ± 15 keV. 

A close connection may exist between the production 
of non-thermal hard X-rays and radio emission. How- 
ever, until recently all AXPs were radio quiet. For AXP 
IE 2259+58 the detection of radio em ission has now 
been claimed by iMalofeev et al ] (120051 ). but the obser- 
vations and analysis were di fficult and confirmation is 
required. lHalpern et all ( 2005| ) discovered transient radio 
emission from the transient AXP XTE J1810-045 which 
appeared sharply modulat ed at the rotation p eriod with 
peak flux densities >1 Jy ( Camilo et al.l l2006) . orders of 
magnitude brighter tha n the reported upper limits for 
this or any other AXP. iGaensler et al.l (|200lh have ob- 
served 4U 0142+61 with the VLA (1.4 GHz), but only a 



2006d) . 

We present the first quasi-simultaneous multi- wavelength 
observation campaign to study 4U 0142+61 from the ra- 
dio up to hard X-rays. 



2 Multi-wavelength campaign 

For June- July 2005, 1 Ms INTEGRAL (20-300 keV) 
dedicated 4U 0142+61 observations were scheduled. With 
these observations, we tried to get, nearly simultane- 
ously, as much wavelength-band coverage of 4U 0142+61 
■as possible. An approved 12 hour radio observation with 
the WSRT (21 cm) was rescheduled to overlap with the 
INTEGRAL observations. A regular RXTE (2-250 keV) 
monitoring observation also fell in the INTEGRAL time 
line. For an X-ray imaging observation a Target of Op- 
portunity (ToO) was granted with Swift (0.2-10 keV). 
Finally, optical and NIR observations were requested and 
approved in the Directors' Discretionary Time (DDT) on 
Gemini North. During two nights 4U 0142+61 was im- 
aged in the K s and r' bands. Unfortunately it was not 
possible to schedule the K s and the INTEGRAL obser- 
vations contemporaneous (see Table [J). 



2.1 Hard X-rays: INTEGRAL 

The INTErnational Gamma-Ray Astrop hysics Labora- 
tory (INTEGRAL; IWinkler et all l2003f ) is ESA's cur- 
rently operational hard X-ray/soft gamma-ray space tele- 
scope. For the study of AXPs t he low-energy detect or of 
IBIS, called ISGRI (20-300 keV lLebrun et al.ll2003l) . has 



Table 1 Multi-wavelength campaign measurements 



Obs 


Time span 


Exposure 


WSRT 


July 2, 2005 


12 h 


Gemini K 3 


July 26, 2005 


1125 s 


Gemini r' 


July 13, 2005 


2400 s 


Swift 


July 11 - 12, 2005 


7400 s 


INTEGRAL 


June 29 - July 17, 2005 


868 ks 
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proven itself to be of great importance. The serendipi- 
tous discovery of AXPs in the INTEGRAL energy band 
was a result of the combination of long exposure times 
and the 29° wide FOV of IBIS. 

For this work we have analysed 1 Ms of observations 
of 4U 0142+61 performed between June 29 and July 17, 
2005. The data were screened for solar flares and erratic 
count rates resulting from passes through the Earth's ra- 
diation belts. After screening the net exposure was 868 
ks (Tabled]). The observations consist of 265 pointings 
(Science Windows, ScWs) which can last up to one hour. 
The ScWs have been analysed separat ely with the offi- 
cial I NTEGRAL software OSA 5.1 (see lGoldwurm et alJ 
120031 for IBIS-ISGRI scientific data analysis) in 20 en- 
ergy bands between 20 keV and 300 keV with exponen- 
tial binning. These analyses result in sky images for ev- 
ery ScW in 20 energy bands. The spectrum was built up 
by averaging the count rates from each ScW, weighted 
by the variance. For the conversion into flux values, the 
spectrum was normalized to the known total Crab spec- 
trum (nebula and pulsar) u s ing a curved power law as 
determined by lKuiper et al ] (|2006h . 

4U 0142+61 is detected up to 230 keV with a 3. Oct 
significance in the 150-230 keV energy band. The total 
spectrum shown in Fig. [T] was fitted with a power-law 
model resulting in a photon index r = 0.79 ± 0.10. The 
quality of the fit is good with a reduced chi square Xr — 
1.12 for 16 degrees of freedom. The 20-230 keV flux is 
(17.0 ± 1.4) x 10" 11 erg cur 2 s _1 . 

Our new total spectrum shows that this AXP is now 
detected at even higher energies than reported earlier, 
without an indica tion for a spectral break. The hint for 
a break found by Iden Hartog et all ( 2006al ) is not con- 
firmed in this data set. 



2.2 Soft X-rays 
2.2.1 Swift 



The Swift-XRT (0.2-10 keV; iBurrows et a"Dl2005l ) obser- 
vation was performed on July 11-12, lasting 8500 s. Of 
this observation 7400 s of data were taken in the Photon- 
Counting mode and were analysed with the FTOOLS 
xrtpipeline, version build-14 under HEADAS 6.0 (Hill 
et al. 2005). Photons were extracted from an annular- 
region (3 pixels inner radius, 30 pixels outer radius) in 
order to avoid pile-up contamination. Background spec- 
tra were taken from close-by source-free regions. 

As mentioned in Sect, [IJ AXP spectra in the X-ray 
domain (< 10 keV) can usually be fitted satisfactorily 
with a black-body and a power-law model. When we use 
this canonical model the fit results are: A^h = (1.01 ± 
0.10) x 10 22 cm 2 ; kT = (0.400±0.012) keV; r = 2.7±0.3; 
Xr = 0.97 (dof = 573). These parameters yield an un- 
absorbed flux of (3.9 ± 0.2) x 10 -10 erg cm -2 s _1 in the 
0.7-6.0 keV band, or (7.8 ±0.5) x 10" 11 erg cm" 2 s" 1 in 



the more commonly used 2-10 keV band. Such a model 
fits the measured X-ray spectrum well, but systemati- 
cally overestimates the flux at the lower X-ray energies. 
Furthermore it seems meaningless for extrapolation to 
the NIR window. In particular, the soft power-law com- 
ponent dominates the black-body component for energies 
less than ~1.5 keV, meaning that also the A^h estimate 
is affected and estimated too high. 

Alternatively, we used a double black-body model 
to fit the measured spectrum, yielding an acceptable 
fit (Xr — 0.98; dof = 573) with an excellent fit at the 
lower X-ray energies, but underestimating the higher X- 
ray energies (> 4.5 keV). The two temperatures are 
0.38 ± 0.02 keV and 0.78 ± 0.10 keV and the N H is 
(0.61 ±0.03) x 10 22 cuT 2 . The unabsorbed flux is (2.11 ± 
0.06) x lO" 10 erg cm" 2 s" 1 in the 0.7-6.0 keV band and 
(6.94±0.29) x 10" 11 erg cur 2 s^ 1 in the 2-10 keV band. 
Arguably, this model is again canonical, but it is more 
appropriate for extrapolation to the NIR, which is shown 
in Fig. [TJ 

2.2.2 Rossi-XTE 

4U 0142+61 moni toring data of the PCA (2-60 keV; 
iJahoda et al.lll996l) on board the Rossi X-ray Timin g Ex- 
plorer (RXTE. iBradt. Rothschild & Swanij|l993l) were 
used to create a phase-coherent timing s olution valid dur- 
ing th e multi- wavelength campaign (see lGavriil fc Kaspil 
2002, for a detailed description and use of the AXP- 
monitoring campaign) . This ep hemeris is essential for IN - 
TEGRAL timing analysis (c.f., Iden Hartog et al.l2006bT ). 
In this work it is used for the WSRT observation in or- 
der to reduce the number of trials in the search for a 
(possible) weak pulsed radio signal. Data from RXTE- 
observation IDs 90076 and 91070 were analysed with the 
pulsar-timing software package TEMPCQ- The resulting 
ephemeris is valid between MJD 53251 and MJD 53619, 
with the following characteristics: Epoch MJD 53420.0, 
v = 0.1150929855(7) Hz, v = -2.639(8) x lO^Hzs- 1 
and v = 3(2) x 10~ 23 Hzs- 2 . 



2.3 Optical & NIR: Gemini 

The field of 4U 0142+61 was ima ged in the K s band w ith 
the Near Infrared Imager (NIRI: lHodapp et al.ll2003h on 
Gemini North, Hawaii, in the night of July 26th, 2005. 
NIRI has standard broad-band and narrow-band filters 
covering 1-5/xm. 

The final image was created by subtracting dark frames 
from each science frame, dividing by a flat field derived 
from the images themselves, and then aligning and stack- 
ing all the images. The photometry was measured using 
the PSF-fitting package DAOPhot (|Stetson|[l987h . and 
calibrated against the photometry provided for several 

1 http://pulsar.princeton.edu/tempo 
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field stars in lHulleman. van Kerkwijk fc Kulkarnil(|2004l ). 
The K s magnitude was found to be 19.96(10). 

Optical r'-band images were obtained on the night 
of July 13th, 2005 using the Gemini-North Multi-Object 
Spectrograph fGMOS-N; lHook et al.ll2004D . We subtracted 
the bias and divided by screen flats, before stacking and 
photometering the images. For the calibration, the pho- 
tometry listed in Hulleman, van Kerkwijk & Kulkarni 
(2004) was used, interpolating between the R an d V- 
bands using the relationship in ISmith et al.l (|2002f ). We 
find for 4U 0142+61 r' = 25.42(6), where the statistical 
uncertainty in the measurement and the uncertainty in 
the calibration of the photometry zero-point are similar. 



2.4 Radio: WSRT 

Using the Westerbork Synthesis Radio Telescope (WSRT) 
we have searched for both pulsed and unpulsed radio 
emission. An observation of 12 hour duration was car- 
ried out at a frequency of 1380 MHz (~ 21 cm) with a 
bandwidth of 80 MHz. Using the synthesis data a map 
was made using standard routines in the MIRIAE0 pack- 
age. In the resulting map (rms ~30 /iJy) no source was 
detected at the location of 4U 0142+61 leading to a 
3cr upper limit on its flux of ~90 /iJy. Simultaneously 
we also summed the signals from all 14 dishes of the 
WSRT coherently to form a so-called tied arr ay. The data 
were then sent to the Pulsar Machine PuMa ( Voute et al.l 
12002ft which formed a digital filter bank with 512 chan- 
nels and a sampling time of 409.6 /is. As the dispersion 
measure in the direction of the source was unknown, we 
tried many trial dispersion measures and then folded 
each one with the RXTE ephemeris (see Sect. I2.2.2p . 
Each resultant profile was then inspected to determine 
if there was a significant detection. We also performed a 
standard pulsar search analysis on the full data set. Nei- 
ther the folding nor the search revealed any significant 
detection of radio pulsations. A 5<r upper limit was de- 
termined at 77 x yj~[d/l — d)/iJy where d is the duty cycle 
of the pulsar. Using d = 0.5, a typical value in the X-ray 
regime, the 3cr upper limit is ~46 /iJy. It has to be noted 
that for this analysis the whole o bservation wa s used . 
A fi ner analysis li k e perf ormed by lHalpern et al (2005) 
and ICamilo et al.l (|2006h who discovered the transient 
AXP XTE J1810-197 as a bright transient radio source 
in smaller intervals is still ongoing. However 4U 0142+61 
has not exhibited the same sort of transient behaviour 
in X-rays like XTE J181 0-197, specially it has not shown 
any large outburst, and therefore the radio characteris- 
tics of these sources might be very different. 



http: / / www.atnf.csiro.au / computing/ software/miriad / 



3 Multi-wavelength Source-Energy Distribution 

In Fig.[T]the quasi-simultaneous multi- wavelength Source- 
Energy Distribution (SED) is presented covering roughly 
10 orders of magnitude in photon energy. In the lower 
left corner the WSRT (1380 MHz) continuum (upper ar- 
row) and timing (lower arrow) 2a upper limits are shown. 
The Gemini K s (2.15 /im) and r' (0.6 /mi) data points 
are dereddened assuming Ay = 3.6 as determined by 
iDurant fc van Kerkwijk! |2006bf) . The Swift-XRT spec- 
trum between 0.7 keV and 6 keV is corrected for absorp- 
tion with ]Vh = 0.61 x 10 22 cm -2 . Also shown is the corre- 
sponding double black-body fit (solid line), extrapolated 
towards lower energies (dotted line). The INTEGRAL 
flux values between 20 keV and 230 keV show the hard 
X-ray spectrum for this AXP. The power-law fit (solid 
line) is also extrapolated to lower energies (dashed line) . 

Significant flux variability has been reported in the 
optica l, NIR and soft X-ray bands ( Durant fc van Kerkwijk! 
2006c). therefore in this quasi-simultaneous spectrum we 
can investigate better how the fluxes in the different 
bands relate to each other. 



4 Discussion 

The results for the different wave-band measurements 
render separately no surprises. The Swift-XRT spectrum 
shows a typical soft AXP spectrum and NIR and opti- 
cal magnitudes are around the earlier reported values. 
The INTEGRAL spectrum is in agreement with previ- 
ously found results, however, thanks to the high effective 
exposure in this dedicated observation the maximum en- 
ergy up to which 4U 0142+61 could be detected is now 
higher, namely 230 keV. The spectrum, described with a 
power-law model with photon index r — 0.79 ± 0.10 and 
luminosity 8.7 x lO^ei-gs^ 1 (20-100 keV, d = 3.6 kpc) 
can be compared with th e spectr al results reported ear- 
lier by Iden Hartog et all (|2006al ) using an independent 



data set: T = 0.73+0.17 an d luminosity 8 5 x 10 i 4 erg s" 1 
(20-100 keV, d = 3.6 kpcV iKuiper et al.l (|2006l ) derived 
r = 1.05 + 0.11 and luminosity 8.1 xlO 34 ergs -1 (20-100 
keV, d — 3.6 kpc). All these findings are within errors in 
agreement. Therefore there is no evidence for long-term 
ti me variability at hard X-ra y energies yet. 

Iden Hartog et all ( 2006af ) published 2a flux upper lim- 
its for 4U 0142+61 in the 0.75-30 MeV window analysing 
COMPTEL data collected over the years 1991-2000. These 
observations are obviously not contemporaneous to the 
multi-wavelengths campaign reported here, but the ap- 
parent stability at hard X-ray energies seems to justify a 
direct comparison of the COMPTEL upper limits with 
the hard X-ray spectrum measured with INTEGRAL. 
Thus assuming that the hard X-ray emission is stable, 
Fig. [T] clearly shows that the hard X-ray / soft gamma- 
ray spectrum of 4U0142+61 has to break between ~ 200 
keV and 750 keV in order not to be in conflict with the 
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Fig. 1 Source-energy distribution of 4U 0142+61 from radio up to gamma-ray energies. For observation details, see Sect. [3] 
The WSRT radio limits are shown for the continuum emission (top arrow) and the pulsed emission (lower arrow) . Extrapo- 
lations to lower energies are shown for the Swift-XRT double black-body fit (dotted line) and for the INTEGRAL power-law 
fit (dashed line). It can be seen that neither the soft X-ray nor the hard X-ray spectral fits extrapolate to the Optical and 
NlR fluxes. An optical-NIR excess orders of magnitude above these extrapolations is evident. The COMPTEL upper limits 
at MeV energies do not belong to the multi-wavelength campaign. 







COMPTEL upper limits. If we were to assume that by 
chance 4U0142+61 was in a low state at MeV energies 
during the long COMPTEL observations, it would be re- 
markable that also for the other AXPs detected similarly 
by INTEGRAL (IE 1841-045 and 1RXS J1708-4009) no 
signal was found at MeV energies during the many ob- 
servations that they w ere in the COMPTE L field-of-view 
spread over 10 years (jKuiper et alj 12006). We conclude 
that a drastic break is required. 

The Swift-XRT and INTEGRAL spectra can not be 
fitted simultaneously with a two-component model con- 
sisting of i.e. a black body plus only one power law com- 
ponent. The excess of high-energy photons (w.r.t a single 
black-body spectrum) in the Swift spectrum can not be 
accounted for by the upcoming hard X-ray spectrum. 
This is evident in Fig. [TJ The extrapolation of the hard 
X-ray spectral fit towards lower energies is already five 
times lower at 6 keV than measured with Swift. 

The relation of the optical and NIR flux values to 
the soft and hard X-ray spectra is still an enigma. It was 
noted earlier (e.g., iHulleman. van Kerkwijk fc Kulkarnil 
l2004h that extrapolations of fits to the soft X-ray spectra 
to the optical and NIR bands are not consistent with the 
measured variable flux values. This is more than evident 
in Fig.[TJ where we show the extrapolation of the double 
black-body fit, which reaches optical and NIR fluxes ~4- 
5 orders of magnitude lower than the measured values. 
The very hard X-ray spectra above 10 keV, originally in- 
terpreted as being of non-thermal origin, led to sugges- 
tions of a common origin as the likely non-thermal opti- 



cal and NIR emissions. However, Fig. [T] unambiguously 
shows that also extrapolation of the power-law spectral 
shape measured at hard X-ray energies passes 2-3 or- 
ders of magnitude underneath the optical and NIR data 
points. The hard X-ray, soft X-ray, optical and NIR com- 
ponents seem to have different origins. 

The discovery of the hard X-ray spectra of AXPs, 
stimulated new attempts to search for non-thermal radio 
emission possibly originating from the same sites and/or 
production mechanisms (e.g. synchrotron radiation) in 
the magnetar's magnetosphere. The WSRT upper limits 
at 1380 MHz are however not constraining in this re- 
spect. The INTEGRAL power-law fit extrapolates also 
orders of magnitudes below the radio upper limits. 

The physical interpretation of the hard X-ray spec- 
trum measured by INTEGRAL is not evident. Assuming 
a distance of 3.6 kpc as recently determined by Durant 
& van Kerkwijk (2006a), the 20-230 keV flux (Sect. 12!]) 
translates to a luminosity of 2.6 x 10 35 erg s _1 , which ex- 
ceeds the maximum luminosity available from rotational 
energy 1OSS0 with a factor of ~2000. Moreover, the hard 
X-ray luminosity is comparable with the high soft X-ray 
luminosity (1.1 x 10 ergs , 2-10 keV) as measured 
with Swift (see Sect. |2~2~T]) . 

Over the last two years theoretical attempts have 
been made trying to explain the new hard X-ray emis- 
sion above 10 keV from AXPs. A particularly interesting 



3 £re = 1.3 x 10 32 ergs 1 assuming a neutron star with 
R = 10 km, M = 1.4M 
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interpretation of the high-energy emission is the appli- 
cation of a magnetar-corona model by Thompson & Be- 
loborodov (2005) and Beloborodov & Thompson (2006a, b 
in this volume). In the last papers, a theoretical model 
is developed explaining the formation of a hot corona 
above the surface of a magnetar, so far the only model 
attempting to explain the multi-wavelength emission. In 
essence, the twisted magnetosphere acts as an accelerator 
that converts the toroidal-field energy to particle kinetic 
energy. They show numerically that the corona self or- 
ganises quickly into a quasi-steady state, with a voltage 
of ~1 GeV along the magnetic field lines. Pair production 
occurs at a rate just enough to feed the electric current. 
Interestingly, the heating rate of the corona is ~10 36 - 
10 37 ergs _1 , in agreement with the observed persistent, 
high-energy output of magnetars. Furthermore, they de- 
duce that the static twist will decay on a timescale of 
1-10 years, setting the scale for time variability to look 
for in the high-energy emission. The transition layer be- 
tween the atmosphere and the corona is the likely source 
of the observed ~100 keV (e.g. bremsstrahlung) emission 
from magnetars. Finally, it is worth mentioning that the 
corona emits curvature radiation which can also supply 
the observed IR-optical luminosity. Of particular con- 
cern in the hot coronae model is the reported timescale 
of 1-10 years. For 4U 0142+61 we do not see evidence 
for strong variability in the hard X-ray emission over 
the first years of INTEGRAL observations, nor over the 
longer period of RXTE monitoring observations. Future 
work will concentrate on this aspect, as well on detailed 
timing studies, including phase resolved spectroscopy. 
Obviously, as much INTEGRAL data as possible will 
be collected on 4U 0142+61 and the other hard X-ray 
emitting AXPs to search for the break energy in the total 
spectrum, a key parameter in the discrimination between 
the different proposed models. More detailed analysis of 
this campaign and additional INTEGRAL o bservations 
will b e pr esented in forthcoming p apers by iRea et ail 
(2006) and lden Hartog et alJ (|2006rJ l. respectively 
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